######################################################
# 
# --- Command log file ---
# PLAXIS 2D 2025.1
# Tutorial - 10
# Excavation and dewatering
# Note: Use the PLAXIS Commands runner to execute this
# 
######################################################
_initializerectangular SoilContour 0 0 100 30
_setproperties "Title" "Excavation and dewatering" "Company" "Bentley Systems Inc" "Comments" "" "UnitForce" "kN" "UnitLength" "m" "UnitTime" "day" "UnitTemperature" "K" "UnitEnergy" "kJ" "UnitPower" "kW" "WaterWeight" 10 "ReferenceTemperature" 293.15 "LiquidSpecificHeatCapacity" 4181.3 "LiquidThermalConductivity" 0.0006 "LiquidLatentHeat" 334000 "LiquidThermalExpansion" 0.00021 "LiquidTemperature" 293.15 "IceSpecificHeatCapacity" 2108 "IceThermalConductivity" 0.00222 "IceThermalExpansion" 5E-5 "VapourSpecificHeatCapacity" 1930 "VapourThermalConductivity" 2.5E-5 "VapourSpecificGasConstant" 461.5 "UseTemperatureDepWaterPropsTable" False "ModelType" "Plane strain" "ElementType" "15-Noded"
_borehole 0
_soillayer 0
_soillayer 0
_soillayer 0
_set Soillayer_1.Zones[0].Top 30
_set Soillayer_1.Zones[0].Bottom 27
_set Soillayer_2.Zones[0].Bottom 15
_set Borehole_1.Head 23
_soilmat "Identification" "Silt" "SoilModel" "Hardening Soil" "Colour" 15262369  "gammaUnsat" 16 "gammaSat" 20 "E50Ref" 20000 "EOedRef" 20000 "EURRef" 60000 "nuUR" 0.2 "PowerM" 0.5 "cRef" 1  "phi" 30 "UseDefaults" True "GroundwaterClassificationType" "USDA"  "GroundwaterSoilClassUSDA" "Silt" "GwUseDefaults" True "GwDefaultsMethod" "From data set" "InterfaceStrengthDetermination" "Manual" "Rinter" 0.65 "GapClosure" True
_soilmat "Identification" "Sand" "SoilModel" "Hardening Soil" "Colour" 10676870 "gammaUnsat" 17 "gammaSat" 20 "E50Ref" 30000 "EOedRef" 30000 "EURRef" 90000 "nuUR" 0.2 "PowerM" 0.5 "phi" 34 "psi" 4 "UseDefaults" True "GroundwaterClassificationType" "USDA" "GroundwaterSoilClassUSDA" "Sand" "GwUseDefaults" True "GwDefaultsMethod" "From data set" "InterfaceStrengthDetermination" "Manual" "Rinter" 0.7 "GapClosure" True
_soilmat "Identification" "Loam" "SoilModel" "Hardening Soil" "Colour" 10283244  "gammaUnsat" 17 "gammaSat" 19 "E50Ref" 12000 "EOedRef" 8000 "EURRef" 36000 "nuUR" 0.2 "PowerM" 0.8 "cRef" 5 "phi" 29 "UseDefaults" True "GroundwaterClassificationType" "USDA" "GroundwaterSoilClassUSDA" "Loam" "GwUseDefaults" True "GwDefaultsMethod" "From data set" "GapClosure" True "POP" 25 
_setmaterial Soillayer_1 Silt
_setmaterial Soillayer_2 Sand
_setmaterial Soillayer_3 Loam
_gotostructures
_plate (40 30) (40 14)
_plate (60 30) (60 14)
_platemat "Identification" "Diaphragm wall" "MaterialType" "Elastic" "Colour" 16711680 "W" 8.3 "PreventPunching" True "EA1" 12000000 "EI" 120000 "StructNu" 0.15 
_setmaterial Plates DiaphragmWall
_posinterface (Line_1 Line_2)
_neginterface (Line_1 Line_2)
_line (40 23) (60 23)
_line (40 20) (60 20)
_n2nanchor (40 27) (31 21)
_n2nanchor (60 27) (69 21)
_n2nanchor Point_5 (31 17)
_n2nanchor Point_6 (69 17)
anchormat "Identification" "Anchor rod" "MaterialType" "Elastic" "EA" 500000  "LSpacing" 2.5
_setmaterial NodeToNodeAnchors AnchorRod
_embeddedbeam Point_10 (28 19)
_embeddedbeam Point_12 (72 19)
_embeddedbeam Point_13 (28 15)
_embeddedbeam Point_14 (72 15)
_embeddedbeammat "Identification" "Grout body" "MaterialType" "Elastic" "Colour" 9392839 "LSpacing" 2.5 "Diameter" 0.3 "E" 7070000 "TSkinStartMax" 400 "TSkinEndMax" 400   
_setmaterial embeddedbeams GroutBody
_set (embeddedbeam_1.Connection embeddedbeam_2.Connection embeddedbeam_3.Connection embeddedbeam_4.Connection) "Free"
_group (Line_5 Line_9 Line_6 Line_10)
_group (Line_7 Line_11 Line_8 Line_12)
_rename Group_1 "GroundAnchor_Top"
_rename Group_2 "GroundAnchor_Bottom"
_lineload (28 30) (38 30)
_well 42 20 42 17
_well 58 20 58 17
_mergeequivalents Geometry
_gotomesh
_set BoreholePolygon_2_5.CoarsenessFactor 0.25
_mesh 0.06
_gotostages
_phase InitialPhase
_set Model.CurrentPhase Phase_1
_set Phase_1.Identification "Activation wall and load"
_activate Plates Phase_1
_activate Interfaces Phase_1
_activate LineLoads Phase_1
_set LineLoad_1_1.qy_start Phase_1 -10
_phase Phase_1
_set Model.CurrentPhase Phase_2
_set Phase_2.Identification "First excavation"
_deactivate (BoreholePolygon_1_2) Phase_2
_phase Phase_2
_set Model.CurrentPhase Phase_3
_set Phase_3.Identification "First anchor row"
_activate GroundAnchor_top Phase_3
_set (NodeToNodeAnchor_1_1.AdjustPrestress NodeToNodeAnchor_2_1.AdjustPrestress) Phase_3 True
_set (NodeToNodeAnchor_1_1.PrestressForce NodeToNodeAnchor_2_1.PrestressForce) Phase_3 500
_phase Phase_3
_set Model.CurrentPhase Phase_4
_set Phase_4.Identification "Second excavation"
_deactivate (BoreholePolygon_2_1) Phase_4
_phase Phase_4
_set Model.CurrentPhase Phase_5
_set Phase_5.Identification "Second anchor row"
_activate GroundAnchor_bottom Phase_5
_set (NodeToNodeAnchor_3_1.AdjustPrestress NodeToNodeAnchor_4_1.AdjustPrestress) Phase_5 True
_set (NodeToNodeAnchor_3_1.PrestressForce NodeToNodeAnchor_4_1.PrestressForce) Phase_5 1000
_phase Phase_5
###########################################################################
#New Dewatering Phase - Modifying Phase 6_Final excavation (Tutorial 6 ) 
###########################################################################
_set Model.CurrentPhase Phase_6
_set Phase_6.Identification "Final excavation"
_set Phase_6.PorePresCalcType "Steady state groundwater flow"
_deactivate (BoreholePolygon_2_2) Phase_6
_gotowater
_waterlevel (0 23) (40 20) (60 20) (100 23) 
_rename UserWaterLevel_1 "LoweredWaterLevel"
_set Water.GlobalWaterLevel Phase_6 LoweredWaterLevel
_gotostages
_activate (Well_1_1 Well_2_1) Phase_6
_set (Well_1_1.Q Well_2_1.Q) Phase_6 1.5
_set (Well_1_1.Hmin Well_2_1.Hmin) Phase_6 17
_selectmeshpoints
/output addcurvepoint "Stresspoint" Soil_2_3 37 18
/output addcurvepoint "Stresspoint" Soil_2_5 46 18
/output update